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We investigate reliability of Gamow- Teller transition strengths computed in the proton-neutron 
random phase approximation, comparing with exact results from diagonalization in full Ghu) shell- 
model spaces. By allowing the Hartree-Fock state to be deformed, we obtain good results for a wide 
variety of nuclides, even though we do not project onto good angular momentum. We suggest that 
deformation is as important or more so than pairing for Gamow- Teller transitions. 



PACS numbers: 21.60. Jz, 21. 60.Cs,23.40.-s 



I. INTRODUCTION 



Weak processes, such as beta (/?) and double-beta 
'(3) decay, have deep consequences for nucleosynthesis 
Ij and physics beyond the Standard Model Q. Weak 
processes are also sensitive to details of nuclear struc- 
ture: allowed Gamow- Teller (GT) transitions depend in 
particular upon PauH blocking 0,0]. Thus, predictions 
for astrophysics as well as interpretation of /3/9-decay ex- 
periments must be modeled with care. Because of this 
sensitivity to Pauli blocking and thus upon configura- 
tion mixing, the large-basis interacting shell model (SM) 
^l] provides one of the best microscopic approaches to 
Gamow- Teller transitions. 



The interacting shell model, however, is computation- 
ally expensive, and only recently has the full QTioj pf- 
shell become tractable. A simpler approach is the ran- 
dom phase approximation (RPA) and its generalizations, 
which have been widely and successfully applied to giant 
resonances such as the electric dipole (giant dipole reso- 
nance or GDR) and has also been applied to many 
important problems in weak transitions 1, 6, 7, 8, 9, lOj. 
It is not immediately obvious, however, that the RPA is 
an appropriate approximation for all transitions. Be- 
cause the RPA is the small-amplitude limit of time- 
dependent mean- field theory prL Il2j | , it seems appropri- 
ate for the GDR, which can be described semiclassically 
in the Goldhaber- Teller model ^3 as protons oscillat- 
ing in bulk against neutrons. The application of the 
RPA and its extensions to Gamow- Teller transitions, al- 
though with a long history more problematic. 
Two implicit assumptions in the RPA are, first, ground- 
state correlations on top of a mean-field state are small, 
and, second, particle-hole phonons have boson commu- 
tation rules, which means that Pauli blocking is not fully 
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treated. Therefore, as GT transitions are sensitive to 
Pauli blocking, they may not be well matched to RPA 
calculations. 

Reading the literature only furthers these doubts. A 
number of authors have previously tested the efficacy 
of calculating GT transitions within the proton-neutron 
quasiparticle RPA (pnQRPA) , through comparison with 
exact calculations either with full shell-model diagonal- 
ization 0,0,0,01 or group-theoretical schematic mod- 
els The efhcacy of these pnQRPA calculations, 
which we will discuss in more detail below, can be broadly 
summarized as poor, typically overestimating the total 
transition strength and underestimating the first moment 
of the transition strength. In order to compare with the 
calculations described below, is important to understand 
that these calculations used spherical J = 0, iV-even Z- 
even parent states and treated pairing carefully, starting 
with either the Hartree-Fock-Bogoliubov or the Bardeen- 
Cooper-Schrieffer equations. 

With the exception of those tests of Gamow- Teller 
transitions, however, most other tests of the RPA and 
its generalizations have been against toy models. Re- 
cently we began to systematically test the RPA against 
full shell-model diagonalization ,19, 2(1. 21]. In this paper 
we describe the generalization to proton-neutron RPA 
(pnRPA) and compare Gamow- Teller transitions against 
the full shell-model results. In contrast to previous ap- 
proaches, we do not treat pairing carefully, but do allow 
arbitrary deformation, even though we do not project 
to good angular momentum. Furthermore we are not 
limited to even-even nuclides. We obtain good transition 
strengths, and where we can compare to published spher- 
ical pnQRPA strengths, our calculations are generally 
superior. We therefore conclude that proper treatment 
of deformation, even without projection, is at least as 
important as proper treatment of pairing, and arguably 
more so. This agrees with shell model studies that show 
a correlation between E2 transition strengths (a measure 
of deformation) and GT transition strengths 0, . 

We also briefly argue that the pnRPA frequencies 
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should be real, and not complex; such an argument is 
missing from the literature. Included in this discussion 
is a lemma helpful to understanding the solution of the 
pnRPA eigenvalue equations. 



II. FORMALISM 

In this paper we consider only change-changing 
Gamow- Teller transitions; the transition operators are 
thus: 



0± = gA^T±, 



(1) 



where r+ changes a neutron into a proton. Because here 
we consider only strength distributions and not absolute 
transition strengths, we drop the axial vector coupling 
gA- The transition strength from the parent state \i) 
(here always a ground state) to a final state / at excita- 
tion energy E = Ef — Ei is given by 



SiE) 



E,-E)\{f\0\^)\' 



(2) 



The transition strength to a given level is also often called 
the B(GT) value. 

A convenient way to characterize the distribution of 
transition strength is through moments. The zeroth mo- 
ment is the total transition strength, 5*0, 



(3) 



E 



while the centroid of the distribution E is the first mo- 
ment 



e = J2es{e)/So 



(4) 



(note that this centroid is relative to the parent energy 
Ei) and the width of the distribution, AE, is the square 
root of the second moment, calculated relative to the 
centroid: 



AE^ = J2{E-EfS{E)/So. 



(5) 



An important check of any calculation is the well- 
known Ikeda sum rule, which holds true for any parent 
state, both in the SM and in the RPA: 



So{r)-Soir)^HN -Z) 



(6) 



We now briefly review the random phase approxima- 
tion and its variant, the proton- neutron RPA (pnRPA) 
[rHI |. Unlike the standard or like-particle RPA, pnRPA 
provides means to calculate excited states of neighboring 



isobars. The starting point, however, is similar to the 
regular RPA: a mean-field solution of the parent nucleus, 
which in turn defines a particle-hole basis. In our case 
we began with a Hartree-Fock state |HF), allowing unre- 
stricted deformation. The RPA can be derived a number 
of ways, but it essentially approximates the energy sur- 
face about the mean-field state as a harmonic oscillator. 
This leads to an RPA ground state |RPA), which implic- 
itly includes zero-point fluctuations about the mean-fleld 
state |HF), and excited states are one-phonon excitations 
on the ground state: 



/3l|RPA). 



(7) 



For RPA one usually assumes the creation operators are 
simple one-particle, one-hole operators. In order to prop- 
erly discuss the pnRPA, it is useful to further, if briefly, 
recapitulate the like-particle RPA In the matrix 

formalism, one solves 



A B 
B A 



(8) 



Here the matrix A, a Hermitian matrix, is in fact the sub- 
matrix of H between one-particle, one-hole states (also 
the matrix for the Tamm-Dancoff approximation), while 
the elements of B, a symmetric matrix, are the matrix 
elements of H between the HF state and two-particle, 
two-hole states. (For simplicity we assume real HF wave- 
functions so that all subsequent quantities are real.) One 
can show that the excitation frequencies f2 are real if the 
stability matrix, 



A B 
B A 



(9) 



has no negative eigenvalues; this is equivalent to the 
Hartree-Fock state being at a (local) minimum. Part 
of the proof of this requires that both A ± B have no 
negative eigenvalues which in turn requires that A 
have no negative eigenvalues. We emphasize this point 
because the situation will be different for pnRPA. 

One can show that the solutions of Eq. come in 
pairs: if {X,Y) is a solution with frequency Q, then 
{Y,X) is also a solution with frequency —il. The so- 
lution with > is associated with the vector {X, Y) 
such that \X\ > \Y\, and one chooses a normalization 
|A"p — |Fp = 1; that one can do this also derives from the 
nonnegative eigenvalues of the stability matrix. The spe- 
cial case 51 = corresponds to invariance of the Hamil- 
tonian under a symmetry, for example, rotation or trans- 
lation; here the vector {X, Y) cannot be normalized, as 
\X\ = \ Y\, and one must resort to a different formalism 
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The proton- neutron RPA is similar to the hke-particle RPA but with important differences. For hke-particle 
RPA, the phonon creation operator uses one-particle, one-hole operators of the form tt^tt, v'^v (thus the name 
like-particle). The pnRPA phonon creation operator is also one-body but changes the third component of isospin T^: 



mi 



The first and the fourth terms in this equation describe the excited states of the (Z-f- 1, A^— 1) isobar, while the second 
and third the (Z — 1, A^ -I- 1) isobar. We follow the standard convention and use indices m, n for 'particle' states, that 
is, unfilled single-particle states in the Hartree-Fock wavefunction, and i,j for 'hole' states, or filled single-particle 
states; thus ir'^i'i destroys a filled neutron state (or creates a neutron hole in the HF wavefunction, in the usual view) 
and creates a proton in an excited particle state. The equation of motion 

(RPA|[[(5/3, [H,f3^]]\RPA) = f7(RPA|[(5/3, /jl'JIRPA), (11) 

transforms as usual to a non-hermitian eigenvalue problem |l2j 




X{pn) \ 
X{np) 
Y{np) 
Y{pn) J 



= n 



I X{pn) \ 
X{np) 
Y{np) 

\ Y{pn) ) 



(12) 



where the definitions for ^p"'"p and _B"p^p" matrices are similar to the regular proton-neutron conserving formalism, 
where one approximates |RPA) « |HF): 

<".:"n" = (HF|[i.^„„ [i?,7rti^,]]|HF) = « - er)5™n% - (13) 



= -(HF|[^t^7r„ [<^„i/]]|HF) = -VZ,^. (14) 

(For the quasiparticle RPA (QRPA) one uses instead of the Hartee-Fock state a Hartree-Fock-Bogoliubov state or 
a Bardeen-Cooper-Schrieffer state, and the QRPA phonon is composed of quasiparticle-quasihole operators instead.) 
The matrices AP'^'"p and BP'^'"p are defined similarly, but are distinct unless Z = A^; in fact, they have different 
dimensions unless Z = N. Let , NJ^ be number of proton particle and hole states, respectively, and N^, the 
number of neutron particle and hole states. Thus the vectors X{pn) and Y{np) are of length N^N^ while vectors 
X{np), Y{pn) are of length N^NJ^; the two lengths are unequal unless Z = N. Similarly, A"^'^" is a square matrix 
of dimension NpNj^ while A^^'^p is a square matrix of dimension NpNJ^, while B"^'''" is a rectangular matrix of 

dimension N^Nj^ x N^NJl, and BP"'"p = (B"P'P")'^ . 



The zeroes in Eq. (|12|) occur because the original 
Hamiltonian conserves charge (and thus Tz). The overall 
form of Eq. H12|) is identical to that of Eq. © ; we have 
merely introduced a block structure to A, B. Because of 
the zeroes in Eq. (|12() . the equations decouple: 




_^pn,np _^pn,np j ^ y {pn) ) ^ \ Y {pjl) ) 

(15) 



_Qnp,pn _^np,pn J y {np) ) ' \ Y {np) J 



(16) 

It is important to note two things here. First, the 
decoupled equations lfT31l and l(TB)l are not of the same 
form as Eq. ©, because AP"'"p ^ a^p-p" unless Z = N. 
Because of this, some of the usual theorems do not imme- 
diately apply to (|15|) and (|16|l . especially regarding the 



positivity of Q (more about this, however, in a moment) 
and the sign of \X\^ - \Y\^. 

Second, any solution (X{pn),Y{pn)) to Eq. (|15|l with 
frequency O is related to a solution {X{np),Y{np)) 
of Eq. itTCI) with frequency — fi, by {X{np),Y{np)) = 
(Y {pn) , X (pn)) , up to some overall phase factor. This 
simplifies finding solutions. Let {U, V) be a solution of 
Eq. IjlSI) with frequency uj, which can be positive or nega- 
tive. If|L/|^ — > 0, and so normalizable to 1, then let 
{X{pn),Y{pn)) — (U,V) with frequency — uj. Other- 
wise, if |C/|2 - |T/|2 < 0, then let {X{np),Y{np)) ^ {V, U) 
with frequency — —uj be the solution to Eq. H16|l . In 
both cases one can normalize to 1. 

What about — 0? In like-particle RPA, if exact 
symmetries are broken by the mean-field, one obtains 
zero eigenvalues; the corresponding eigenvectors, which 
identify the generators of the broken symmetries [23j, 
are not normalizable and must be treated with care. For 
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pnRPA, although zero eigenvalues are not excluded, the 
corresponding eigenvectors do not play a special role and 
are in fact normalizable, as argued below. 

In like-particle RPA, one can show that stability of 
the mean-field state implies that > 0. But now in 
pnRPA one can have < 0, even for \X\'^ — \Y\'^ = 1. 
At first glance this is troubling; however, we provide two 
arguments which resolve this apparent paradox. 

First, because pnRPA allows charge-changing phonons 
of the form (|HJ|) . one should, at least implicitly, also per- 
form the Hartree-Fock minimization allowing mixing of 
proton and neutron states, with N-Z fixed by a Lagrange 
multiplier A. (We have in fact done such a calculation, 
but the results are indistinguishable from fixing N-Z by 
hand, that is, as one varies A, N-Z makes integral jumps. 
Nonetheless, the image is useful.) Now the constrained 
Hartree-Fock is at a true minimum, and if one considers 
the eigenfrequencies of (|12|l all the resultant pnRPA fre- 
quencies are proven positive-definite (although in prac- 
tice one solves (|15|l or (|16ll instead). To interpret the 
results as physical transition energies, however, one has 
to subtract off A(A^ — Z) which leads to negative frequen- 
cies. 

One can see this directly. By adding the Langrange 
multiplier, one shifts the diagonals of the A matrices, 

i"f 'P" = 'P" - A(iV - Z) • /, (17) 

and 

= + \{N -Z)-I. (18) 

By substituting into Eq. (|15|l . it is clear one is adding a 
constant along the entire diagonal, which only shifts to 
~ O + A(A^ — Z), while the frequencies in Eq. H16|) are 
shifted in the opposite direction. In fact, for any value of 



A one only shifts the frequencies and the solutions (X, Y) 
are manifestly unchanged. 

We call this result the shift lemma: by adding the 
Lagrange multiplier one shifts the relative position of 
the proton and neutron Fermi surfaces, but the only re- 
sult is shifting the absolute value of the pnRPA frequen- 
cies. The relative frequencies and the eigenvectors are 
unchanged. There is a useful consequence: if one obtains 
a zero mode, from the shift lemma one can find a case 
where the frequency corresponding to the same eigenvec- 
tor is positive, and one expects the pnRPA vector to be 
normalizable; and as the eigenvector is independent of 
the shift, it is always always normalizable. The reader 
should note that the shift lemma only arises because of 
the unique isospin dependence of the block decomposi- 
tion (|12|l : no such shift is possible for like-particle RPA 
because one cannot add a constant to the diagonal. 



III. RESULTS 

We test the reliability of pnRPA's predictions for GT 
strengths against exact diagonalization in full QHuj SM 
spaces. That is, we calculate the GT strength distri- 
butions for several nuclei in the sd (lsi/2-0d3/2-0c?5/2) 
shell on top of an inert ^^O core, and two Ti isotopes 
in the pf (IP1/2-IP3/2-O/5/2-O/7/2) sheh above a ^"Ca 
core. While we do not compare directly with experiment, 
we use in our calculations phenomenological interactions 
which are very successful in reproducing the experimental 
data: Wildenthal 'USD' in the sd shell jH], and Richter- 
Brown in the pf shell . The shell model diagonaliza- 
tion were performed using a descendent of the Glasgow 
code [23, and the shell- model strength distributions us- 
ing an efficient Lanczos moment method |28|. 



Having defined the pn-phonon creation operator (3'^ in 110(1 , one can calculate the transition matrix element required 
by Eq. © as 

{RVk\0\\z±iMTi)) = (RPA|[0,/3l]|RPA) = {Xi,{pn/np)0,m + Yy{pn/np)0,^) . (19) 

mi 

In our case, O is the GT transition operator defined in Q|, which induces transitions between the correlated ground 
state of the parent nucleus and the (Z -I- 1, iV — 1) or {Z — 1 , + 1) isobars. 



Table summarizes results for total transition 
strengths, centroids and distribution widths, where in- 
deed we find that the pnRPA moments are reasonably 
close to SM. As a check, note that the Ikeda sum rule 
© is fulfilled in both SM and pnRPA, as expected. We 
find same features typical of like-particle RPA, that is the 
pnRPA centroids are usually lower in energy, while the 
SM distribution widths are larger |^ . The latter can 
be understood as particle-hole correlation beyond RPA 
which further fragment the distributions. 



Note that we have results not only for even-Z, even-N 
nuclides but also odd-odd and odd-A, all with compara- 
ble success. (A technical note: In our mean-field calcu- 
lations, we allow only real wave functions. This does not 
have any effect for the even-even nuclei. For even-odd or 
odd-A nuclei however, the exact mean-field solution could 
be complex, and we see small symptoms of this restric- 
tion: because the rotations with respect to x and z axes 
are complex, the proton and neutron number conserving 
RPA formalism identifies some generators of the broken 
symmetries as lying at small, but not zero, excitation en- 
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TABLE I: Total strength So, centroid E, and width A£ for 
GT transition operator. The nuclei have been grouped into 
even-even, odd-odd and oAA-A. 



So E (MeV) AE (MeV) 



Nucleus 






SM 


RPA 


SM 


RPA 


SM 


RPA 


2°Ne 






0.55 


0.69 


15.81 


12.20 


4.22 


2.42 


22Ne 






0.50 


0.63 


19.71 


16.17 


3.81 


1.33 








6.50 


6.63 


4.48 


4.75 


5.64 


3.79 


2*Ne 






0.51 


0.61 


19.73 


18.37 


3.31 


1.36 








12.51 


12.61 


3.82 


4.26 


4.91 


3.46 


2^Mg 


/? 




2.33 


2.73 


13.40 


10.92 


3.86 


2.33 


2«Mg 






1.78 


2.05 


15.72 


13.42 


3.55 


1.97 








7.78 


8.05 


6.94 


5.93 


5.58 


4.62 


28Si 


/? 




3.89 


3.39 


13.54 


12.29 


3.07 


1.86 


30Si 






2.52 


2.33 


15.59 


13.39 


2.59 


1.83 








8.52 


8.33 


8.69 


7.38 


4.69 


3.26 


32g 






4.01 


4.25 


12.48 


10.38 


3.04 


2.21 


34g 






1.59 


1.88 


14.22 


11.90 


2.53 


2.13 






/?- 


7.59 


7.88 


7.91 


7.88 


4.06 


2.34 


36 Ar 




+//3- 


2.10 


2.22 


12.09 


10.07 


2.59 


2.57 




P 


V/3- 


0.61 


0.79 


9.95 


7.96 


2.27 


1.50 


46^i 




/3+ 


0.44 


0.60 


12.47 


10.46 


1.80 


0.55 








6.44 


6.60 


2.99 


3.31 


3.51 


2.39 


24Na 




/3+ 


1.67 


1.92 


14.59 


12.34 


3.53 


2.65 






/3- 


7.67 


7.92 


6.67 


6.15 


4.87 


3.72 


26^1 


P 




4.28 


4.28 


11.86 


10.37 


3.43 


2.87 


2iNe 






0.63 


0.67 


15.85 


13.96 


4.49 


2.82 








3.63 


3.67 


6.49 


5.82 


5.05 


4.17 


^^Na 






1.39 


1.50 


15.96 


14.06 


3.27 


1.95 








10.39 


10.50 


5.27 


4.92 


5.13 


3.97 


2^A1 






3.20 


2.52 


14.04 


12.72 


2.94 


2.02 








6.20 


5.52 


9.33 


8.61 


5.20 


3.75 


28 Al 






1.80 


1.97 


16.10 


13.51 


2.76 


2.04 






/3" 


10.80 


10.97 


6.61 


6.45 


4.85 


2.99 



ergies. While these approximations can be relevant for 
other transitions '23|, they do not have an impact here.) 

Figs. illustrate selected results in detail. These 
are typical results, neither better nor worse on the av- 
erage. These figures follow a useful convention used by 
many many authors and plot the accumulated sum of 
the strength, Y.B{GT) = J2E<EfS{E) , which allows 
one to compare by eye the first few moments. Again we 
see by eye generally good results, for even-even, odd-odd, 
and odd-A alike. The major systematic error of RPA ap- 
pears to be lower centroids for transitions; why there 
is no similar lowering of the P~ centroid is unknown. 

In Figs.nand|21wc compare /3~^ calculation against the 
spherical pnQRPA calculations of Lauritzen 16], which 
are similar to those of 0, 0| . (All of these papers com- 
puted only (3~^ decay, and only for even-even nuclides.) 
In general the QRPA strengths So are about twice as 



much as the exact calculation, and the centroids E are 
significantly lower. The calculations of found that 
pnQRPA gives results very similar to shell-model calcu- 
lations restricted to 2 particle, 2 hole excitation out of 
the 0(1^/2 orbit. This is believeable, as RPA correlations 
are approximately 2 particle, 2 hole in nature. By con- 
trast, our pnRPA calculations on top of a deformed HF 
state do much better. In particular, our deformed pn- 
RPA calculations better approximate the correct total 
strength than does the spherical pnQRPA. That defor- 
mation can "quench" Gamow- Teller strength is already 
known from shell model calculations 0,122. (See, how- 
ever, pnQRPA calculations of j^^, who find the total 
strength to be insensitive to deformation, although the 
details of the distribution they find to be sensitive. They 
speculate that the difference is that they include up to 10 
harmonic oscillator shells in their calculations, while our 
calculations and those of I3| and l25| are within a single 
harmonic oscillator shell. Why a multi-shell calculation 
should be insensitive to deformation is not clear.) 

There is a curious exception: the spherical pnQRPA 
calculation of 0| on ^^Mg yields So and E with ap- 
proximately the same accuracy relative to SM results 
as us. They only compute the one f3~^ decay, how- 
ever, and it is difficult to understand the difference be- 
tween their methodology and that of Hi^^, 18]; further- 
more, although in principle all three calculations are us- 
ing the same shell- model interaction , which we prop- 
erly scaled to A = 26, we can reproduce Lauritzen's shell 
model results 0] but not those of [T^ . 

In order to further illustrate how important it is to 
treat correctly deformation, we show in Fig. [S] three 
different pnRPA calculations for "^^S: with dashed line 
we represent the distribution strength for a spherical HF 
state, with dotted line, a distribution obtained by start- 
ing with an oblate HF state and with dash-dotted line a 
distribution obtained starting with a triaxial HF state. 
The two deformed states are almost degenerate in en- 
ergy, the triaxial state being slightly lower. From Fig. [S] 
one sees that the strength distribution off the triaxial HF 
state best approximates the SM result (as one might ex- 
pect). QRPA strength distributions which use a starting 
spherical mean field solution IT^ are slightly more frag- 
mented than in our spherical calculation (we obtain just 
one state with J = 1), but the total strength is about 
the same in both pnRPA and QRPA, overestimating se- 
riously the SM total strength. We conclude that using 
a deformed mean-field solution is at least as important, 
and arguably more so, than treating pairing with rigor. 

In the above discussion we have focussed on the gross 
properties of the Gamow- Teller "resonance," which is the 
main application of RPA. But for cold systems, one is 
often most interested exclusively in low-lying transitions. 
As is observable in the figures, our pnRPA calculations 
are rather mediocre when it comes to the lowest-lying GT 
transitions: arguably better than pnQRPA, but not very 
good compared to the exact shell model results. This is 
an important caveat for any applications. 
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FIG. 1: SM (full curve), RPA (dashed), and QRPA(dotted) summed strength for 13+ (left) GT transitions from ^^Ne, ^"S. The 
QRPA calculation is from fl^ 
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FIG. 2: SM (full curve) and RPA (dashed curve) summed strength for (3+ (left) and /3 (right) GT transitions in ^'^Mg. For 
P+ transitions we also include the QRPA calculation [l^ (dotted curve) 




E (MeV) 

FIG. 3: SM (full curve) and RPA (dashed curve) summed 
IV. CONCLUSIONS 

The main purpose of this paper was to investigate 
the pnRPA's reliability for predicting GT transition 
strengths; the motivation is astrophysical applications, 
where, aside from binding energies and electomagnetic 
transitions, a good knowledge of week transitions is es- 
sential. Our tool in this investigation was the interacting 
SM which can provide the exact numerical solution in a 
restricted space. Although we made our tests for nuclides 
near the bottom of the valley of stability, presumably 
our results apply out to the driplines; in fact the major 
uncertainty will be the shell-model interaction, not the 
pnRPA. 




E (MeV) 



strength for /3+ (left) and (3 (right) GT transitions in ^""Na. 

We show a very good agreement between SM and pn- 
RPA transition distributions in a large number of 
nuclides in the sd and pf shells, similar to high-lying col- 
lective electromagnetic transitions investigated in a pre- 
vious paper 21] . Furthermore, we obtain better results 
than spherical pnQRPA, obtaining better suppression of 
the total strength. This may seem surprising, as we do 
not treat rigorously the pairing interaction; on the other 
hand, we model correctly the deformation and this proves 
suitable for describing GT transition strengths. (Because 
of this, we view tests of pnRPA in schematic models with- 
out symmetry breaking j30| as being inadequate.) It is 
possible that a deformed pnQRPA approach, with better 
treatment of pairing, will lead to further improvement. 
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E (MeV) 



FIG. 5: Comparison of the SM (full curve) and pnRPA predictions for the summed strength of GT transitions in S. The 
pnRPA was performed on top of a spherical (dotted curve), oblate (dash-dotted), and triaxial (dashed) HF state respectively, 
emphasizing the importance of a correct treatment of deformation for a correct description of the strength distribution. 



We hope however to extend our program to the deformed 
pnQRPA in the not-so-distant future. Furthermore, it 
will be important to consider multi-shell spaces, where 
at least one set of deformed pnQRPA calculations (not 
validated by direct shell model calculations, however, due 
to the enormous size of the model space) suggest a weaker 
dependence on deformation (2^. 
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